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Abstract. In this paper we consider the question of observational signatures of a false 
vacuum decay event in the early universe followed by a period of inflation; in particular, 
motivated by the string landscape, we consider decays in which the parent vacuum has a 
smaller number of large dimensions than the current vacuum, which leads to an anisotropic 
universe. We go beyond previous studies, and examine the effects on the CMB temperature 
and polarisation power spectra, due to both scalar and tensor modes, and consider not only 
late-time effects but also the full cosmological perturbation theory at early times. We find 
that whilst the scalar mode behaves as one would expect, and the effects of anisotropy at 
early times are sub-dominant to the late-time effects already studied, for the tensor modes 
in fact the the early-time effects grow with multipole and can become much larger than 
one would expect, even dominating over the late-time effects. Thus these effects should be 
included if one is looking for such a signal in the tensor modes. 


1 Introduction 


The study of false vacuum decay in quantum field theory is a topic with a long history, dating 
back to the initial work of Callan, Coleman, and de-Luccia [1-3]; the basic idea is that if the 
theory has multiple vacuum states, and initially is in one with higher energy, henceforth the 
false vacuum, then it is possible to tunnel to a state in which there is a ‘bubble’ of lower- 
energy (henceforth, true) vacuum sitting in the remaining region of false vacuum, the two 
regions being separated by a transitionary bubble wall. This bubble then proceeds to grow 
in size until, in the asymptotically far future, the whole universe is in the true vacuum. 

Within the bubble it is possible to construct the coordinate system of an open FLRW 
universe (with the bubble wall at infinity) [3], which leads to the exciting possibility that if 
our universe were open, it could have arisen out of such a false vacuum decay event, followed 
by a short period of inflation (in part to dilute away the spatial curvature, which would 
dominate initially). Much work was done concerning this idea in the early ‘90’s, for example 
see [4-7], and there has since been a resurgence of interest due the idea of the multiverse, 
for example [8-12], since if such a landscape of possible metastable vacuum states exists 
(e.g. [13-15]), then tunnelling events between these states would seem inevitable. Thus a 
tunnelling event in the origin of our universe is, in some sense, a prediction of the landscape 
hypothesis, and it behoves us to investigate the observational consequences. 

In particular, this paper examines the consequences of tunnelling events in which the 
number of large dimensions in the false vacuum is smaller than the number of large dimensions 
in the true vacuum, as could occur for example if the tunnelling field is the modulus of one 
(or more) of the dimensions. This possibility is motivated in part by string theory; the string 
landscape arises out of the myriad ways to compactify any of the ten dimensions required 
for the consistency of the theory (including the possible flux quanta that can exist on closed 
sub-manifolds, i.e., homology ‘cycles’ of various dimensionality, of the compactification), in 
order to generate a low energy effective theory in a smaller number of dimensions. Including 
the many discrete choices of fluxes there are generically vastly more ways to compactify a 
larger number of dimensions than a smaller number. This leads one to expect that there 
are more vacua in the landscape which have a smaller number of large, apparently non¬ 
compact, spatial dimensions. (See [16, 17] for attempts to tackle the measure problem in 
such a scenario.) Thus in the the (string) landscape a tunnelling event which increases the 
number of large dimensions is a motivated and intriguing possibility [18-22]. This is not to 
say that such events are exclusive to string theory, for example even the Standard Model 
may support lower dimensional vacua [23, 24]. 

Finally, many other approaches to quantum gravity besides string theory suggest that 
space-time may be of lower dimensionality at high energies [25, 26]. This could manifest 
itself observationally in many ways [27, 28], and in particular in a scenario similar to that 
considered here, with inflation after a dimensional transition [29]. 

Such dimension-changing false vacuum decay events and their observational signatures 
have been considered before. The key point is the that the universe that results will have 
large background anisotropy immediately after its creation (at least in the 2 -|- 1 —?• 3 -|- 1 and 
1 -|-1 —> 3 -|-1 cases). The background anisotropy reasserts itself during the matter dominated 
epoch, and the effect of this on the CMB temperature power spectrum was examined in [30] , 
while in [31, 32], and more recently [33], the tunnelling process itself was investigated, along 
with the effect on the scalar power spectrum in the limit of no metric fluctuations. In this 
paper we extend this to the polarisation power spectra, specifically going beyond previous 
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work on the effects of early-time anisotropy by performing full cosmological perturbation 


theory allowing us to examine the effect on the tensor perturbation modes. 

The paper is structured as follows: section 2 discusses in more detail the model which 
we are considering, as well as the background evolution of the universe; section 3 considers 
perturbation theory in a universe with background anisotropy, derives the primordial power 
spectra for the case of2-|-l —)• 3-1-1, and analyses how to obtain CMB power spectra, 
ignoring the effect of late-time anisotropy, which is then considered in section 4. Section 
5 then gathers together the results, compares the effects, and comments on observability. 
Finally section 6 concludes. Four appendices contain more detail on the various calculations 
presented in the main body of the paper. 

2 An anisotropic universe 

We here focus on the case of2-|-l —)• 3-1-1 and 1-1-1 —)• 3-1-1 tunnelling events resulting in our 
present universe but with anisotropy of a particular form. The universe inside the bubble (our 
universe) is spatially homogenous, but anisotropic, with the possibility of curvature in two 
of the spatial dimensions and cosmic shear between those and the third spatial dimension; 
its metric is 



( 2 . 1 ) 


The sign of k depends on the origin of the {r,4>) dimensions. If they are ‘old’, i.e. the 2 
dimension is the one which has recently decompactified (so the transition was 2-|-l —)• 3-1-1), 


then K = —1, since the coordinate system inside the bubble (with the ‘new’ z dimension 
suppressed) must be open. In this case (2.1) describes a Bianchi III universe. On the other 


hand if (r, 4>) are the new dimensions, and the transition was 1-1-1 —?■ 3-1-1, then k G {—1, 0,1} 


depending on the details of their previous compactification (e.g. compactification on 5^ yields 


K = 1). In this case (2.1) describes a Bianchi III, Bianchi I, or Kantowski-Sachs universe 
respectively. 

From the tt and zz components of the Einstein equations one finds 



( 2 . 2 ) 


(2.3) 


which are then supplemented with the equation of motion for matter, or the rr or (pcf) com¬ 
ponent of the Einstein equations. Here p and pz are the tt and components respectively 
of the stress tensor Mpi is the reduced Planck mass, and = (a/a)^, H'^ = (b/b)'^. 
We shall assume that inflation is driven by a single scalar field (p, with potential V, in which 
case during inflation the equations become 



(2.4) 


(2.5) 


0 + (2H, + Hh)^ + V' = 0. 


( 2 . 6 ) 
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The initial conditions depend on the origin of the dimensions, with the scale factor of the ‘old’ 
dimensions starting from zero, and we assume that the scale factor for the new dimensions 
starts at rest from some finite value (which would be expected if the tunnelling field were the 
modulus for the extra dimension) along with the inflaton. 

The evolution of the background has been studied in [30] and will be briefly recounted 
here. Given the boundary conditions mentioned above, immediately after nucleation of the 
bubble the behaviour depends on which dimensions are new: in the 2 + 1 —?■ 3 + 1 case we 
have a{t) « t + 0{t^) and b{t) ^ 6(0) + C>(^^)^ 5 the 1 + 1 —>■ 3 + 1 case 

depends on the sign of the curvature. As inflation proceeds the anisotropy in expansion rate 
will be driven to zero, as will the anisotropic curvature parameter 

In the radiation and matter dominated epochs, like any curvature, this will begin to grow, 
and will source a corresponding anisotropy in the expansion rate (in [30] it was shown that 

//a ^ 

It is important to note that there are two types of anisotropy in (2.1): shear anisotropy 
due to the fact that Ha 7 ^ and curvature anisotropy due to the fact that only two of the 
three spatial dimensions are curved. Thus the observational consequences of this model can 
arise both from anisotropy at early times, which we examine in section 3, and from anisotropy 
at late times, which we discuss in section 4. 

3 Early-time anisotropy 

The effect of anisotropy at early times is to modify the behaviour of the perturbations to the 
inflaton and metric, whose fluctuations will eventually be imprinted upon the CMB. Thus 
we need to extend the usual cosmological perturbation theory to deal with universes of the 
form (2.1); in order to do this we will closely follow [34] in which this was done for universes 
with anisotropic shear, but no curvature. (See also [35] for the same, and [36] in which is 
examined the case of anisotropic curvature, but no shear.) 

We immediately note two things: 

Since we have two scale factors, it is not obvious how we should define conformal time, 
and in fact the appropriate definition will depend on the specific case (i.e. which are the 
new dimensions) under consideration; thus for now we just define c{t)dr] = dt, where 
c{t) is some function of a and 6 . 

As we no longer have a residual SO{3) symmetry on the spatial slices, we cannot use the 
action of this to classify the perturbations, and instead must use the 50(2) symmetry 
of (r, (/)); thus the perturbations will be 50(2) scalars (1 degree of freedom) and vectors 
(also 1 d.o.f. since it has two components and is required to be transverse).^ 

The gory details of anisotropic perturbation theory are relegated to appendix A, and here 
we simply state the results. Of the 50(2) vectors only one survives, and comparison with 

^Note that SO{2) tensors have four components and are required to be symmetric, traceless, and transverse, 
which eliminates all the d.o.f. and thus they do not appear. 


(ii) 


(i) 


(2.7) 
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the isotropic case^ reveals this to be the analogue of the tensor perturbation , and so we 
denote it likewise. The quadratic action for it is given by 

^ / dr] d^Q dk3y/j{\h'^\‘^. (3.1) 

Here —q^ is the eigenvalue of the 2-dimensional Laplacian^ on (r, 0), — A:| is the eigenvalue 
of the 1-dimensional Laplacian for the 2 ; dimension, 7 is the determinant of the conformal 
(r, 4>) metric, and the effective, conformal-time-dependent mass is given by 



Note that in the isotropic limit this reproduces the expected result, — \- 

Two SO{2) scalar perturbations survive, which are the analogues of the isotropic scalar 
mode V (the Mukhanov-Sasaki variable) and the tensor mode /i+. These two are coupled and 
their action takes the form 

= ^ f dq d^q dks^ (h'^KH' + H'^KH - H^MH^ , (3.4) 


where = ((/?, S), and the matrices K, K, and M are given in Appendix A. In order to find 
the precise relationship between {(p, S}, and {u, /i+}, the mixing in the hrst two terms must 
be removed through a suitable field redefinition. Whilst this can be done in principle, and 
it is possible to do it in closed form in the absence of curvature [34], the correct procedure 
for the general case eludes us here. Fortunately, in the simplihed case we consider below the 
demixing will happen automatically. 

We can now also answer the question of how to define conformal time. The requirement 
is that, for each mode, the adiabatic vacuum state (the generalisation of usual vacuum state 
in Minkowski space). 


u{q) 



■i Lo{‘r]')dri' 

? 


(3.5) 


is a solution of the equations of motion as t —)■ 0. This corresponds to uo' <C and > 0 
as t —)• 0. By inspecting (3.3) we see that this is satisfied when c is the scale factor of the old 
dimensions (i.e. those for which the scale factor goes to zero as t does), and the same is true 
for the modes in (3.4). 


^Note that due to the presence of anisotropic curvature, the isotropic limit is not just a = b, but also 
requires sending fife —>■ 0. 

^Note that it satisfies the following relations depending on the value of k [37]: 


f^2 + | ifft — —1, 
<1^ = \ k2 if K = 0, 

[z(/ -1-1) if AC = 1, 


where A :2 G R and Z G N®. Though for compactness we will also write = fcf when ac = 1. 


(3.2) 
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3.1 Primordial power spectra 

We are interested in the power spectrum of the fluctuations at the end of inflation. Since in¬ 
flation efficiently isotropises the background, the modes become their isotropic counterparts; 
therefore in the remainder of the paper we will refer to them as such. The relation of the 
primordial power spectra to the expectation value of the modes is the same as in the isotropic 
case, i.e. for the comoving curvature perturbation TZ = ^v, 

acj) 

(7^(k)7^*(k')) = ^<5(k-k')P7e(fe,//), (3.6) 

and similarly for the tensor perturbations '0-i-,x = Here + A;| (see footnote 

3 for the definition of /C 2 ) and /r = k • z is the cosine of the angle that k makes with the z 
direction. The fact that the power spectrum only has this dependence on the direction of 
the wavevector is due to the residual isotropy in the (r, 4>) plane. 

3.2 CMB correlators 

The fact that the primordial power spectra are not solely functions of k means that their rela¬ 
tions to the correlators for temperature and polarisation fluctuations of the CMB take more 
complicated expressions. In this section we state and explain the results, but relegate their 
explicit derivation to Appendix B. Furthermore we ignore the effect of late-time anisotropy, 
with which we will deal in section 4. 


3.2.1 Scalar mode 


For the comoving curvature perturbation this has already been calculated in [34], and we 
have 


^{S)XY 



dkAS)x 
k ^ 


{kAr^)A\P^{kArj)Pl,fl{k), 


(3.7) 


where X,Y £ {T,E}, is the transfer function for temperature or polarisation,^ Arj = 

Vo - VLSS, and 




{21 + 1){2V + 1){1 - - m)l 

{I + m)\{l' -\- m)\ 


'-i 


dMH^(/x)P^(/a)H7^(/c,/i). (3.8) 


Note that when P-ji depends only on k, due to the orthogonality of the associated Legendre 
polynomials we have Phi^ik) = 2SiiiP'ji{k), and the usual isotropic result is recovered. 

If Pti is an even function of /x (as will be the case here, since parity is still preserved), 
then only multipoles separated by an even number in Z-space will be correlated. Also note 
that the diagonality in m arrises from choosing our isotropic coordinate system to align with 
the anisotropic one (along with the fact that there is residual isotropy in the (r, </>) plane), 
which in general would not be the case and would require a rotation. 


3.2.2 Tensor modes 

Unsurprisingly the case for the tensor modes is significantly more complicated due to the 
fact that the tensor mode functions have spin weight 2. Here we follow the formalism of 
[38], though the derivation in this case is complicated by the fact that we cannot arbitrarily 

^These are listed in Appendix B, in table 1. 
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rotate our coordinate system, due to the anisotropy of -P+,x- The details again can be found 
in Appendix B, and the result is 


MT)XY 

^Wmm' 



dk (T)X 
k ^ 


{kA7j)A[P^{kArj)Pi;!'l{k), 


(3.9) 


where 


l {2l + l){2V + l)il-m)\ik-m)\ 
Y {I + + m)l 

l—\m—2\) P — \m—2\) 


E E PilmPi'l'm [ (3.10) 

i=0 i'=0 d-l 


and 

P^^\k,p) = {l + axaY{-iy-''+P{P+ + Px) + i-iy{cTx + (TY{-iy-''^P{P+-Px), (3.11) 


where cty^e = 1 and as = —1, and P is the parity of the primordial power spectrum (even 
in this case); the constant coefficients f^um are given in Appendix B. 

We have the sum 


l—\m—2\ 

E ft 

i=0 


2{l — i + m — 2)! 


(2(1 — i) + 1)(Z — i — m + 2)! 


= doAl 


2{l + 2)\{l + m)l 
{21 + 1){1 - 2)1(1 - m)l 


(3.12) 


and hence the usual isotropic result is recovered when P^P depends only on k. 

For modes of the same parity (e.g. TE) we again find that only modes separated by an 
even number in Z-space are correlated. However, intriguingly, for modes of opposite parity 
(e.g. TB), for which there is no correlation in the isotropic case, here there is the possibility 
for correlation of modes separated by an odd number in Z-space. Such a signal is interesting 
because it is solely due to early-time anisotropy (whereas the other signals discussed also 
receive contributions from late-time anisotropy). 


3.3 Quasi-de Sitter approximation 

In order to solve the equations of motion for the perturbation mode functions, and thus find 
the primordial power spectra, we would have to resort to numerical methods in the general 
case, however they simplify dramatically if we take a quasi-de Sitter approximation in which 
the energy density is only due to vacuum energy (i.e. the inflaton is frozen). In fact the 
modes in (3.4) decouple, and for the canonically normalised scalar mode we are left with 


LO 


2 

V 


= C 




(3.13) 


The expression for is significantly more complicated, and so we do not show it here, but 
instead it can be found in Appendix A. In the isotropic limit these both reduce to Zc^ — 
as expected for mode functions on de-Sitter space. 
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For the rest of this section we will specialise to the case in which the z dimension is 
the new one (so k = — 1 and c = a) - performing the analysis for the other cases should be 
straightforward. In this case solving the equations of motion for the scale factors gives 


a{t) = 


' 3Mpi w ^ 

^ cosech(—??), b[t) = 


—^coth(-r/), 


(3.14) 


where the conformal time is r] = Intanh ^2 ' Plugging this into (3.13) yields 


= A:^(l — tanh^ ^) “ 2 cosech^ ^ ~ 4 


(3.15) 


and through an appropriate set of field redefinitions the equation resulting from this can 
be recast into the hypergeometric differential equation (see the Appendix C for the details). 
Thus we can solve it, and evaluate the result at the end of inflation, 77 —)■ 0, to get 


P-Jl oc 


27r^ 


lim — 

5 ?-s-o a 


sinh^ irk + cos^ vr/c/x 


Gdvr^ sinh irk 

1 cosh Trk — sin nk^ 


1 1 


F —- + -/c/x + -ik r —7 — + -ik 


1 1 


Stt^ 


sinh Tik 


l + lk-^2^?-l) + 0[k-^)], 


(3.16) 

(3.17) 

(3.18) 


where k = k\/l — /x^ = A;sin0. This is plotted as a function of k, for various values of 
/X = COS0 in figure 1. We see that it oscillates as a function of k when the wavevector has 
a component along the z direction; mathematically this is due to the cos( 7 rfc 7 x) term in the 
numerator of (3.17), and physically it is due to the fact that the comoving particle horizon 
in the z direction is finite (since it starts with finite size), which gives rise to resonance 
effects. The magnitude grows as the wavevector approaches the new z dimension, due to the 
sinh( 7 rA:) term in the denominator of (3.17), and ultimately to the fact that the adiabatic 
vacuum initial condition blows up there. 

The exact effect on the CMB correlators will be discussed in section 5, but for now let 
us note that from (3.18) we see that for large k the behaviour can be roughly split into two 
parts, depending on the value oi k = k sin 0 (the projection of the wavevector along the old 
dimensions): 


Pn 


Po(l + |A: ^cos^ 0) for ksm.6>kt, 

Po^-^ for ksine<h, 


(3.19) 


for some transition value kt ^ below which Pti rises sharply. This rise can be attributed to 
the adiabatic vacuum initial condition we have used, which also blows up in this regime. In 
reality, immediately after tunnelling the mode function will be in a more complicated state 
than the adiabatic vacuum, see for example [12, 39], and thus the small k rise in the power 
spectrum should perhaps be taken with a pinch of salt. In section 5 we shall see that most 
of the interesting signatures are not in fact too dependent on this. 

Now we turn to the hx mode. Plugging (3.14) into (3.3) yields 


2 __2 , \ + k\l+A^?) 

^ ^ (I + fc 2 (l —/x^ tanh^ r/)) cosh^ 77 


3A;^7x^(| + A;^) 

[(1 + ^^(1 “ fP tanh^r/)) cosh^ 77 ]^ 


(3.20) 
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Figure 1. The primordial power spectrum for the coming curvature perturbation as a function of 
wavevector {k = 1 corresponds to the curvature scale), for various values of the angle 9 between the 
wavevector and new 2 dimension. The peak around k = 1 followed by oscillatory behaviour when 9 is 
small, is due to the fact that the comoving particle horizon in the 0 direction is finite (since it starts 
with finite size), which gives rise to resonance effects. 


The extra terms relative to the scalar mode are only important when k < 1, i.e. when the 
wavevector is pointing significantly in the direction of the new dimension, and so outside this 
regime we should be safe using the first line of (3.19). The same is true for the /i_|_ mode, 
though we do not show here due to its length. Outside of this regime we will need to use 
the full forms in (3.20) and (A. 16), however numerical studies we have done indicate that 
the behaviour of the tensor power spectra at small k is similar to that of Pn-^ 

The singularity structure of the equation arising from (3.20) is the same as that of Heun’s 
differential equation and thus it is transformable into that, and in appendix C we discuss 
this and attempt to find the power spectrum Px ■ Unfortunately, due to the complexity of 
Heun functions, we are not able to find a complete, closed form expression, however for small 
values of k it does show the k~^ behaviour exhibited in the scalar mode case, encouraging 
us to take their behaviour to be the same as (3.19) (though with the transition value kt as a 
free parameter). 

4 Late-time anisotropy 

The background also acquires a non-negligible anisotropy at late times that grows during 
the periods of radiation and matter domination, and thus we should also include effects due 
to this when studying the impact on the CMB. In [30] this was studied extensively for the 
temperature perturbations, and so here we simply quote their findings. Three effects were 
found: 

(i) the surface of last scattering (LSS) is warped, being ellipsoidal in shape rather than 
spherical; 

(ii) the angle on the LSS from which a given photon originates is offset from the angle at 
which it is received; 

®More detailed numerical studies are required to properly investigate the departures of Px,+ from Pn and 
each other. 
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(iii) the redshift to the LSS is now angle dependent. 


The effect of the last two is to shift the coefficients of the multipoles, so that 

1 


0>lm ^ 0,1m ^kO ^ ^ 2i,m! 


(4.1) 


i=—l 


for some coefficients h\ which are explained in appendix D. The first has the effect that the 
distance to the surface of last scattering acquires an angle dependence: 

X = A:A?7/cAry ^1 - , (4.2) 

and thus the argument of the transfer functions, /cAry, does too. One deals with this (to first 
order in which is small) by Taylor expanding: 

Ai{x) A;(xo) - A((xo)xo^Y^flfco42o(6')- (4.3) 

The derivatives of the transfer functions can be found in appendix D, in table 1. 

We now need to extend this to deal with polarisation, which is particularly easy, since 
the same effects apply in that case apart from the angle-dependent redshift, and thus the 
coefficients h in (4.1) will be different for polarisation than temperature.® 

Combing all the effects we have for both scalar and tensor modes, to first order in 


^XY _ Omm‘ 


dk 

TT Jo k 
1 


AfAj;pi 


ll'm 


+ ^ (^45 V ^^^+2i ~ K ’ ^ J) ) Pl,l'-2i,m 


(4.4) 


where the coefficients / arise from the Y 20 spherical harmonic in (4.3) and can be found in 
appendix D. This is appropriate both when early time anisotropy effects are not included, 
and when they are, with the difference being the expression for Pu'm (i-e. (3.8) and (3.10), 
or their isotropic simplihcations). 


5 Results 

Having examined the various ways in which the background anisotropy due to the tunnelling 
event will affect the CMB correlators, let us now look at the size of the effects and consider 
their observability. We specialise to the case ofa2-|-l—)-3-|-l transition, and use the 
quasi-de Sitter approximation discussed in section 3.3. 

®Note that this effect will still appear in temperature-polarisation cross-correlation, but only through the 
temperature transfer functions. 
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5.1 Observables 


The key observational difference compared to the standard cosmological model is that certain 
off-diagonal elements of the correlation matrix are non-zero, and thus we consider the ratio 
of these to the diagonal elements, 


Cll'm 


and 


Cu'm _ Cw 
'Enroll m Cu' 


(5.1) 


The reason to include the latter is that our results implicitly rely on the coordinate system 
used to calculate the CMB correlators being aligned with the direction of anisotropy, which 
in practice is unknown (though of course one could rotate one’s coordinate system to find an 
estimate of the direction of anisotropy). 

The leading order effects are all proportional to n^o, the value of the curvature parameter 
today. This can be directly constrained in the same way as isotropic curvature, however there 
is an even stronger indirect constraint [30]. The fact that in the TT spectrum modes with 
Al = 2 are correlated means that the CMB temperature quadrupole receives a contribution 
from the monopole of order ToQ^o] since the temperature fluctuations are ~ 10^ times smaller 
than the monopole, must be correspondingly small in order that we do not observe an 
anomalously high quadrupole.' 


5.2 Scalar mode 

Let us first briefly estimate the size of the effect. From (3.19) we find 

OO \ 

^ “ sinTrfc) -F 

i=—oo / 

(5.2) 

for some coefficients c, whose precise values aren’t needed in the following paragraph. To find 
the CMB correlators we then need to integrate over the transfer functions which all consist 
of linear combinations of spherical Bessel functions (possibly multiplied by inverse powers of 
kAr]), and therefore we need to use the integral® dkk~''ji{kAr])jii{kAr]) oc (Ar/)^, and so 
the effects due to early-time anisotropy should be ~ (A?])^. Noting that Ary 2y/VLko/ zlss-, 
we see that the early-time anisotropy effects should be ~ 10“®r2fco- This should be compared 
to the late time anisotropy for which we expect effects ~ flfco when \l — l'\ =2. 

Figure 2 shows the above mentioned observables in units of for each of the two 
classes of effects (early- and late-time anisotropy). We see that both effects are smaller 
than one would expect, but that the late-time effect grows with I, whilst the early-time effect 
decreases. Thus it would seem that, as expected, for the scalar mode for Al = 2 the late-time 
effects completely dominate (as predicted in [30]). 

There are, however, observables for which the early-time effects could come dominate: 
Al = 2n for n > 1. The late-time contribution comes from keeping terms higher order in 
in the expansions in section 4, and so we would expect them to be ~ including higher 

order terms in (3.19) will generate the desired effect, but at the cost of larger inverse powers 
of k which results in larger powers of once the k integral is performed, and thus the 

^This doesn’t, however, give us a precise bound since there can always be cancellations between the different 
contributions to the quadrupole. 

®Since the 1 — sinTrfc factor oscillates much more rapidly than the Bessel functions, we can just take its 
average value. 


Pn,ll'm OC (5//,/ -h k 


-2 
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Figure 2. Al = 2 elements of CMB multipole correlation matrix (compared to the AZ = 0 elements, 
in units of the anisotropic curvature parameter) due to the scalar mode, comparing the late and 
early-time effects; top: m = 0; bottom: sum over all m. The early time effects are sub-dominant, and 
in fact a couple of orders of magnitude smaller than expected 


early-time contribution when the wavevector is not aligned along the new dimension [k > kt) 
will behave likewise. On the other hand the early-time contribution when the wavevector 
is aligned along the new dimension [k < kt) will generate non-zero correlations for Al > 2 
which we see from (5.2) will not come at the expense of higher powers of flfco- 

Looking at the actual coefficient in front of these terms once the k integral is performed, 
and we find that 

^ IQ-^nkon-"^ (5.3) 

Oz,/,o 

(for the I which maximises the effect), which is the desired non-exponential falloff with 
n, and we find that this is similar to or larger than the late-time anisotropy effects for 
n > 2. Unfortunately, the effects we are considering are now extremely small ( 10 ® for 

Qko = 10“^). As we shall see in the next section, the tensor modes may offer a better 
opportunity to study the effects of early-time anisotropy. 
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Figure 3. Al = 2 elements of BB multipole correlation matrix (compared to the = 0 elements, in 
units of the anisotropic curvature parameter), comparing the late and early-time effects; left: m = ±2; 
right: sum over all m. For the early time effects the tensor power spectra are taken to be (3.19), with 
kt, the value of fc = fcsind at which the power spectrum experiences a sharp rise, a free parameter; 
kt = 0 corresponds to no sharp rise, kt = \ to the same behaviour Pji, and fc* = 1 to the earliest 
possible onset of the rise. The early time effects grow quickly with I becoming several orders of 
magnitude larger than than expected, and potentially dominant over the late-time effects. 


5.3 Tensor modes 

As explained in section 3.3, for the two tensor modes although we have been unable to derive 
a complete solution for the primordial power spectra in the quasi-de Sitter approximation, 
when the wavevector is not aligned along the new dimension (k > kt) their behaviour is the 
same as the scalar mode, and there are reasons to believe that their behaviour is similar even 
outside of this regime. Therefore in this section, when examining the early-time anisotropy 
effects we take the two tensor mode power spectra to be equal, and to have the same form 
as the scalar mode power spectrum (but with kt as a free parameter). 

Figure 3 compares for the BB correlators, the late-time anisotropy effects with the early 
time anisotropy effect for various values of kt- It shows a stark difference to the scalar mode. 
We see that, somewhat regardless of the value of kt the early-time effect grows with /, and 
does so faster than the late-time effect, becoming dominant for Z ~ 10 — 40,® even despite the 
factor of -2^55 which bedevils the early-time effect. That this behaviour is present even when 
kt = 0, i.e. when there is no sharp rise in the power spectrum when the wavevector is aligned 
along the new dimension, indicates that this is not simply due to the initial conditions we 
have taken for the modes but is a genuine effect (see the discussion towards the end of section 
3.3). 

This behaviour is also seen in the tensor TT correlators; whilst interesting, this is 
observationally inaccessible due to the the dominance of the scalar mode. For the TE, ET, 
and EE correlators the late-time anisotropy effect grows more quickly (the same is true for 
the scalar mode: see figure 2) and so always dominates over the early-time effect. 

®We stop the graph at Z = 30 because above this factors beyond simply the Sachs-Wolfe effect (which is all 
we have included) come into play when calculating CMB correlators. 
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Figure 4. Al = 1 elements oi TB and EB multipole correlation matrix (compared to the geometric 
mean of the Al = 0 (tensor) TT/EE and BB elements, in units of the anisotropic curvature param¬ 
eter); left: m = ±2; right: sum over all m. The tensor power spectra are taken to be (3.19), with 
kt, the value of fc = fcsin0 at which the power spectrum experiences a sharp rise, a free parameter; 
kt = 0 corresponds to no sharp rise, fc* = j to the same behaviour Pji, and fct = 1 to the earliest 
possible onset of the rise. 


Now we come to a further observational signature, which was previously mentioned in 
section 3.2.2: the possibility oi TB and EB correlations. Parity remains unbroken however, 
and so such correlations can only occur between multipoles for which Al is odd. To investigate 
the size of this effect we again compare to a AZ = 0 correlator in the isotropic case, and as 

the TB and EB correlators are zero, we choose to compare to , where X 

is T or E. It should also be noted that this effect is purely due to early time anisotropy. 

The results can be found in figure 4, and we first note that there is dependence on kt, i.e. 
on the behaviour of the primordial tensor power spectra when the wavevector is aligned along 
the new dimension; if the power spectrum does not experience a sharp rise {kt = 0) then the 
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= 1 EB and TB correlators grow quite slowly, with reaching 

~ 10“^Ofco for I = 30, which is about an order of magnitude smaller than the CPj^ 2 /^a^ 
effect. However if fet > 1/4 (the value for the scalar mode), then the EB and TB effects can 
be competitive with the BB, Al = 2 signal. Unfortunately without having a precise form for 
the tensor mode primordial power spectra we are unable to determine this effect precisely 
and judge how interesting it is. 

6 Conclusions 

In this paper we have examined some of the observational consequences for CMB correlation 
functions, of a model in which our universe arises out of a tunnelling event from a parent 
vacuum with a smaller number of large dimensions. Such a scenario is motivated by the 
string landscape, in which one might expect there to be more vacua with a larger number of 
compactihed (and hence smaller number of large) dimensions. 

We have gone beyond previous work [30-33] by considering not only the temperature, 
but also the polarisation power spectra, and also by performing full cosmological perturbation 
theory at early times in order to investigate the effect not only on the scalar, but also the 
tensor perturbation modes. 

When analysing the size of the various effects we have specialised to the case of a 
2 + 1 —)• 3 + 1 dimension changing pattern; doing the same for the 1 + 1 —)• 3 + 1 case 
should be simple, however it would still be interesting to check if the results are qualitatively 
the same in that case. In the limit in which the energy density of the universe only comes 
from vacuum energy and (anisotropic) spatial curvature we have solved exactly for the mode 
function that becomes the scalar mode at late times; this is also the same as for those which 
become the tensor modes, whenever the wavevector of the perturbation is not aligned along 
the ‘new’ dimension. 

For the scalar mode this has allowed us to confirm what was expected, that the effects 
of early-time anisotropy are sub-dominant to the late-time effects (~ lO^^flfco versus ~ 
lO^^rifco)) fo^^ polarisation as well as temperature, and in fact we have shown that the effects 
are even smaller by a few orders of magnitude than one would naively expect. 

Intriguingly for the tensor modes, however, we have found that for the temperature and 
H-mode spectra the early-time effects grow with multipole, such that they can become equal 
to the late-time effects at around I ~ 30; this behaviour is independent of the precise form 
which we assume for the primordial power spectra when the wavevector is aligned along the 
‘new’ dimension. It would be very interesting to extend this analysis to larger I and see to 
what extent the effects continue to grow, especially their size around I ~ 100. We have shown 
as well that this model produces TB and EB correlations (for multipoles separated by an 
odd value), although the signal in this case is rather small ( 10 ^rifco)- 

All of this shows that not only is it testable, a model in which our universe arises out 
of tunnelling from a vacuum with a smaller number of large dimensions, but also that the 
effects of the anisotropic primordial power spectrum can also be observed in principle. Un¬ 
fortunately, given the current and predicted future limits on the tensor-to-scalar ratio [40] it 
seems unlikely that such a signal would be observable in practice. Given the future prospects 
for observations of the large-scale distribution of matter in the universe with experiements 
such as LSST [41], Euclid [42], WFIRST [43], and SKA [44], and especially Hydrogen 21-cm 
line intensity mapping [45, 46], it would be very interesting to examine what the effects of 
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the model would be on large-scale structure (though again it seems likely that the effects of 
the anisotropic primordial power spectrum would be sub-dominant). Indirectly these things 
thus allow us to probe the landscape hypothesis. 

The author is especially grateful to John March-Russell for suggesting the idea behind this 
work, and for helpful comments on a draft, as well as to Pedro Ferreira for comments. JHCS 
is funded by the STFC. 


A Anisotropic perturbation theory 


In this appendix we cover the details of perturbation theory in a universe with background 
anisotropy both in expansion rate and curvature, but with a residual isotropy in two of the 
spatial directions. The discussion closely follows [34] in which this was done for the case of 
anisotropic expansion rates, but no spatial curvature. 

We assume that inflation is driven by a single, minimally-coupled scalar field, and 
expand around the background: (j) = ^ + (p. Imposing the background field 

equations and neglecting boundary terms, the action at quadratic order in the perturbations 
is 






(V)' - 2^' ) 


+ J 




- (Va V)2) 

(A.l) 


Standard procedure would then be to classify the contributions to hfj,y according to their 
transformation properties under 50(3), but as mentioned in section 3, due to the symmetries 
of this problem we instead classify them according to transformation under 50(2). Similarly 
we leave flexible the definition of conformal time, so that the background metric takes the 
form gfj_y = diag (—c^, 6^) , where 'yij (i,j = 1,2) is the metric on the two dimensional 

space of residual isotropy, whose covariant derivative will be denoted by simply Vj. The 
metric perturbation is decomposed as 

c^ViB + c^Bi \ 

hf^y = j (—2'I'7jj-|-2V(jVj)i?-|-2V(jii'j)) id^B + Cyd^Bi 1 , (A.2) 

V -622S / 

where (4', T, S, B, B, E, x) are scalar perturbations, and (Rj, Bi, Ei) are vector perturbations 
(with two components each); the vectors are required to be transverse (ViB^ = 0 etc.), which 
means that they each carry one degree of freedom, and hence we have ten degrees of freedom 
in total. The derivatives with respect to the third coordinate are chosen for later convenience, 
and the factors in front of the off-diagonal elements, cs etc., are left free for now (though we 
shall see later that they drop out of the final expressions). 
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Let us consider how gauge transformations act on the perturbations. Considering —)> 

and decomposing d^^), we have 


c 

r2 /,2 

^X 


Ei Ei 


^ ^ ^ + -f, 

a 

r2 n2 


Bi y Bi „ 

_4 

ip- (p'f. 


s^s + ^e° + 534, 

B . B - —f 

Cy 

E^E-C, 


(A.3) 


From this we can see which combinations of perturbations can be set to zero, through ap¬ 
propriate gauge transformations. For example, the gauge can be completely fixed by setting 
E = Ei = B = "^ = 0, as will be done in the following. 


A.l Vector modes 

As usual the modes with different transformation properties will not mix at quadratic order, 
and we first focus on the vector modes, of which Bi and Bi remain after fixing the gauge. 
The action (neglecting boundary terms) is 





Bi,3Bi^3 


+ 2c)/ 





2 


a[ir\ 

a b ) 


H—2 


^CyBi^3 (A - k) Bi^3 - ^CyBiABi 



(A.4) 


where A is the two dimensional Laplacian on (r, (/>). We expand in eigenfunctions of this 
(denoting the eigenvalue — and Fourier transform in the z dimension and notice that Bi 
appears without time derivatives; its equation of motion gives 


_ 

WJ %q^ + kl 


kl 




Substituting this into (A.4), and making the field redefinition 


h 


X 


Mpi Cy 





Bi, 


(A.5) 


(A.6) 


then yields (3.1) after some algebra. 

^°This choice is partially motivated by a desire to avoid the curved two dimensional subspace. 
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A.2 Scalar modes 

Of the scalar modes, S, B, x, and ip remain after fixing the gauge; their action is 


+ 4 (J [clB - clxY 




^P 1 L 

i+i+v“\p 


62 


P2 P2 
Cs Q 


(? 0 ? 


y 


y 




« 4 9^ d 2 


(? a?' 


2 a' . kl a' 


-4 ^S + ^-S'-2^^-xhi> 

\ a 0 ^ a 


(A.7) 


X, and B only appear as auxiliary fields, and upon integrating them out one reaches (3.4), 
in which the matrix appearing in the kinetic term is 


a?b ( 2k/3 


1 - 


2re 


\2M3 


0 2KHa<p 

S q-^5^ 


where fd — —V + ^ — 4i?2-|^j and 5 — Ha + Hh + 2H, 
in (3.4) is 


02 , 2M2,//2 

52- H 


4«^M2,H2 

^252 


(A. 8 ) 


The second matrix appearing 


/ 


A = 


2a?b 


1 + 


2k/3 


-1 


M 2 


4> 


koO? \ 




-4k A, ^ 


+ kV' 


-iP{qH^+qHHa-Hb)5 

+K 


-21^32^ /3 + ( 3 A, + y')6 


g 2 M 2 j 


^2 + 2/3 


K5-2HaY, 


-2Kf{2Ha + F') - 20(^ + A,6) 


-2(A,-Ab)(/3 + A,6) 

(A.9) 

Due to the complexity of the effective mass matrix, M, we do not show it here. 

The kinetic part of this can be diagonalised (so as to properly identify the mode func¬ 
tions) by performing the field redefinition h = CH, where the matrix C satisfies the following 
conditions 


C^C = K, 


and 


c^c' - c'^c = ]^{k-k^ 


The action is then 


4ci = 2 / drjd'^Qdks^ I h'^h' - k 


C 


-it 


M+-k'] C-^ - C''C-^ 


h 


(A.IO) 


(All) 


There are enough conditions to uniquely determine C (up to an overall constant unitary 
transformation), and in the absence of spatial curvature one can solve this, however in the 
general case it seems that one may need to resort to numerical methods. 
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A.3 Relation to isotropic modes 

In order for the correct interpretation we must relate these modes to their counterparts in 
the isotropic case {a = b, k, = 0). In longitudinal gauge the metric perturbation takes the 
following form (for a perturbation with wavevector (( 7 , 0 , fes)): 




/- 2 $" 


V 




and the normalised mode functions are 


V = a [ ip + 


0 

^3 riso 


0 

_ qks uiso 
k'^ 

_2/,iso 
- 2 ^'- + 




h+ = 


h — ah}^° 

nx — ^ an^ . 


(A. 12 ) 


(A.13) 


From (A.3) we see that (A. 12 ) can be transformed into (A. 2 ) (in the isotropic limit) in our 
gauge (E = Ei = B = ^ = 0) via the following gauge transformation 

^ 3g^ + 2 fc| 

2q^k‘^ + ■ 

(A.14) 

(Note that due to the transverse condition on SO{2) vectors {kiV^ = 0), and our choice of 
wavevector, in this section we have = 0 for any transverse 50(2) vector.) Applying this we 
can relate the dynamical modes to the isotropic modes and hence to the normalised modes: 

.Mpi Qh ^ /A icN 






Ha. = - 


Mpi 

a/ 2 “A: 2 ^’ 


hx — 


^ a' 


\I/i' 




^ qk ^ ’ 


^ = - 


f + ‘^k‘1 

2 g 2 fe 2 ' 




Thus we have the correct physical interpretation of our modes in the isotropic limit (which 
is valid at the end of inflation). In particular notice that the final expression agrees with the 
isotropic limit of (A. 6 ) (up to a phase), and the hrst two to that which one would get from 
(A. 8 ) in the quasi-de Sitter limit used in the main text. 


A.4 /i+ mode function 

For completeness we here present the effective mass term for the mode function in the 
quasi-de Sitter limit: 


OJ 


2 

+ 


= c 



/ K A7i c6 \ q^ 


n _ 





2i3k \ ^ 

5252 J 


aPbiBi 

c^<52 



2/3k \ ^ 

J 


+ X 


- 2AW ( 1 - ^ 


VEi 


cS 

n, 


1 + 


2/3k 


-1 


a^bUl 

c3<52 


1 - % 


1 _j_ 

^ g 252 


-1 


(A.16) 


where AB = Ba — Bb (the calligraphic font represents conformal time Hubble factors), and 
/3 and 5 are defined below (A. 8 ). 
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B CMB correlators with anisotropic primordial power spectra 


The standard formulae for CMB correlators in terms of primordial power spectra make 
use of the isotropy of the spectra and so must be generalised for our case. The following 
two subsections sketch the derivations of (3.7) and (3.9); the transfer functions and their 
derivatives can be found in table 1. 


scalar 


tensor 


T 

E 

T 

E 

B 


A/(x) 




{^3i{x) - ji+i{x)) 

^ {^ji{x) - ji+i{x)) 

3i{x) - (l - ji+i{x) 


ji{x) 


(1-2)1 ^ 

^ji(x 

^ _ (Z+l)(^+2) \ 


2^^J Jl(x) + ^jl+i(x) (^ + l)(y- 

^jl(x) - ji+i(x) 




Table 1. Transfer functions and their derivatives; derivatives acting directly on spherical Bessel 
functions have been removed using the recursion relation (x) = ^ji{x)—ji+i{x) = ji-i{x) —^^ji{x). 


B.l Scalar mode 


For the scalar mode this is not so difficult, and in fact the appropriate expression for TT 
can be found in [34]. The perturbations projected onto the sky are related to the coming 
curvature perturbation TZ by [38] 

oc = dvr ^ i^{2L + l)f^^^jL{x)YlMYLM{n)n{k) (B.l) 

L,M 

A(^)® oc = dvr ^ i^{2L + jL{x)YlMYLM{n)n{k), (B.2) 

L,M 


where x = kArj and the operators are = 1, £^^') = (1 + In the usual way the 

CMB perturbations are decomposed into spherical harmonics with coefficients aim, and we 
have 


MS)xy 


{a 


(S)X iS)Y 
Im ^I'm' 


d^k d^k! 

(27r)^ (27r)^ 


dfIfidr!ft,y;^(h)yy^,(fi')(A(^)^(k,n)A(^)^(k',h')*) 


(B.3) 


oc 


(4.)V-'''(2L+ 1 )( 2 L'+ 1 ) / f 

/T T f A/f' 

d£l^Y{^{n)YLM{n) 


L,M,L',M' 

^ ‘^%Pn{k)YlMYL^M'{i^) 


{2tt 


d£li,,Yvm'{i^)Yl,M'{i^) 

(B.4) 


where and we have used the definition of the power spectrum in (3.6). The 

angular integrals over h, h' can be trivially performed using the orthogonality of the spherical 
harmonics; due to residual isotropy in the (r, cp) plane, the power spectrum is independent 
of the azimuthal angle of k and thus the integral over this variable can easily be performed 
by rewriting the spherical harmonics using Yim{0,4>) = ~ cos0)e*'"‘^, 

yielding a term proportional to 6mm'- Thus, when all is said and done, one arrives at (3.7). 
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B.2 Tensor modes 

The tensor modes are necessarily more complicated as the tensor perturbations have spin- 
weight 2 and hence are not rotationally invariant. We follow the formalism of [38]; the 
perturbations projected onto the sky are 

oc ^(1 - jl^) ((V^+ - + (^+ +iV’x)e-2*^) (B.5) 

oc ^(1 - jx^) -h (V'-r + g(T)^ixk-n 

oc ^(1 - jl^) ((V^+ - #x)e2'‘^ - (^+ +iV’x)e-2*^) (B.7) 

where x = kAr], the operators are = 1, = (—12 -|- x^(l — — 8xdx), = 

(8x -|- 2x‘^dx), and (/2 = cos 6 = k • h, 0) are the coordinates of the direction of observation, 
n, in a coordinate system in which k = z. Then in the usual way we have 


^i:r)XY _ , iT)X IT)Y*. _ 
— \^lm ^I'm' / “ 

1 f dk dVt 


d^k d?\^' 
(27r)^ (27r)3 


(B.8) 


2 


k (27r)" 


{(p++Px) [(^/L+) {y (iPm'+y) + (y 

+ (P+ - Px) [ux [xipm-) {y {lum'+y) + {^1^+) (3> (^ZW-)*)] } 

(B.9) 


where (Tt,_e = 1 and as = —1, and 




2'i ixk-fi ±210 


(B.IO) 


Note that Ip^ _ = (—1)^+™ j , and so in fact we will only work with and drop 

the ± subscript, and we also take m > 0 since Ip_^ = • 

In the isotropic case these angular integrals are simple to perform as we are free to 
rotate k such that it is along the z direction, and the coordinates of n are simply 
However in the absence of isotropy of P+,x we cannot do this and must keep the directions 
of k and n general. Decomposing the plane wave, we have 


/,^ = 47r^i^(2L + l)ji(x)y£M(k) / dn^YC^{h)YLMim - 

L,M 

and a little geometry reveals 


, 2*0 = (1 _ ^ 2^-1 


1 


-1 + 3c„)sfc - 2snSkCn [cfc cos((/>„ - (/>fc) + ism{4>n - (pk)] 


+ Sj- 


1-7:4 + ick sin (2((0„ - cpk)) 


(B.ll) 


(B.12) 
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where c„ = cosOn, Sn = sin0„, with {On, (pn) the coordinates of h, and similarly for k. Plug¬ 
ging this into (B.ll) and rewriting the spherical harmonics in terms of associated Legendre 
polynomials, we have 



t:{21 + !)(/ - m)\ ^(2L+l)i^jz.(x) ^ (^-(^ + P))! ;- , ^ f ( 


{I -|- m)! 


p=-2 


{L + {m + '’(Mfc)e 


where 


(B.13) 


io,i,L,n. = 

I±2,l,L,m = df, (1 - ^,^)Pr{^^)P^^\^^), 


fo = l-fil (B.14) 

/±1 =-2Y'T^(/ifc ^ 1) (B.15) 
f±2 = l{lTf^kf- (B.16) 


Using recurrence relations for the associated Legendre polynomials, and their orthogonality 
properties, the integrals can be evaluated exactly. Doing this and examining the form of 
for various values of I and m reveals that it takes a surprisingly simple form: 

I ~ ' ~ri ^ l—\m—2\ 

where the coefficients fd are given by 


fdoim = l{l - l){l - m + 1){1 - m + 2 ), 

A>o,«m = (-l)*+^ 2 m + 1 )*(^ - ( 2 / + 1 )) + {m-l) + l{l + 1 )^ , 


2 l*l( 2 (Z + i) - l)( 2 (^ + i) + l)( 2 (/ + i)+ 


(B.18) 

(B.19) 


(B.20) 


Finally, plugging (B.17) into (B.2) and writing X {^ji{x)) = Af(x), one recovers (3.9). 


C Solving for the mode functions in the quasi-de Sitter approximation 
C.l Scalar mode 

In order to find the behaviour of the scalar mode function, and hence the power spectrum of 
the comoving curvature perturbation, we need to solve the following equation 

v" + ^A:^(l — /J.^ tanh^ v) ~ cosech^ ^ ~ \ v = 0, (C.l) 

where rj —)• —oo corresponds to the tunnelling event, and r/ —)• 0 to the end of inflation. 
Changing variables to T = tanh^ r] and <h = \/T{1 — T)P^v, where k = k^l — fp, turns this 
into 

r(i - r)ci>" + (^-1 + (^-1 + r) - (-\k^ + 1^1 + ^^^ = o, (c. 2 ) 
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which we recognise as the hypergeometric equation, and thus the general solution to (C.l) is 


v{ri) = (sech^r/) 2 **^ 


Cl ^ 1 1 ~ 1 1 1 ,2 

j^2Fi(-j + 29.-j-2<i;5;tanh , 

, ^ /5 1,5 1 ,, 5 ,2 

+ C2tanh??2i^i ( 4 + 2^’ 4 “ 2^ ’ 2’^^“ ^ 


(C.3) 


where q = k^ + ik. Setting this equal to the adiabatic vacuum at early times then allows us 
to fix the constants and we find 


Cl = 


„ifcin 2 p 1 4 _ ig^ p 1 _ COS (nkfi) + i sinh |^7rA: 
r(-i)rfifc') 2i sinh (irk 


(C.4) 


The value of C 2 is unimportant as we have 

2 


jU3 


lim - 
v^o a 


2'k‘^ 


lim iv sinh ry) 
)?—^0 


27r^ 


Cl 


(C.5) 


since 2 ^ 1 ( 0 , 6 ; c;0) = 1. This gives (3.17), and to get (3.18) the gamma functions must be 
expanded using 


r( 5 ; + a) a-bS^ Gi{a,b) 


i=0 




T{z + b) 

where Go = 1 , Gi = ^{a — b){a + b — 1), and G 2 = ^ (*^ 2 ^) (3(® + b — 1 )^ — (a — 6 + 1 )). 


(C.6) 


C.2 Tensor modes 

As mentioned in section 3.3, when k > 1 the equations for /i+,x take the same form as 
that for V. The full forms for (3.20) and (A. 16) are complicated, yet it turns out that the 
equation for hx can be transformed into Heun’s differential equation. This can be seen by 
first performing the transformations of the previous subsection, and seeing that the resulting 
equation has regular singularities at T G {0,1, a, 00 }, where a = y^(l + 4 ^)) which is 
precisely the singularity structure of Heun’s equation. Performing the further transformation 
h ^ {a — Tyh with c = ^ ^1 + ^ i then explicitly transforms it into this. 

The general solution would then be a sum of two Heun functions with appropriate 
parameters; to apply the boundary conditions at ry = — 00 T = 1, we need to evaluate 

the Heun function at argument unity, which requires a connection formula relating solutions 
regular around T = 0 to those regular around T = 1. Unfortunately no such formula is 
known for the general case and so we cannot provide a full solution for the power spectrum; 
however some results are known [47] which allows us to say that the behaviour for A: <C 1 is 
the same as for the scalar mode, i.e. Px oc k~^. 


D Coefficients and functions appearing in the late-time anisotropy effects 

Here we state the coefficients appearing in (4.4), the expression for the effects of late-time 
anisotropy. For more details see [30]. 
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The warping of the last-scattering surface leads to the (previously isotropic) transfer 
functions taking on terms which are proportional to Y 20 , and since each transfer function is 
multiplied by Yim in the formula for Cu/mm' these lead to the following coefficients: 


k2okzm 


1 


2,m? 


where 

Am _ 1 I {I - m - 1){1 - m){l + m - 1){1 + m) 

~ V nl2l - ly (2/ -3)(2/ + 1) 

Am _ f + 1 - 3nP 

“ V ^2{2l - l){2l + 3) 

Am _ 1 I {I - m + 1){1 - m + 2){l + m + 1){1 + m + 2) 

~ V nl2l + 3\j (2/ + l)(2Z + 5) 
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The offset of the angle on the LSS at which a given photon is received generates terms 
which look like 

1 

sm{2e)deYUe,<i,) = E (D.5) 

i=—l 

where the coefficients are 

= -^(/ + C = -4/r, 5^ = llfr- (D.6) 

If the angle dependence of the redshift is not included, i.e. we are dealing with polari¬ 
sation, then the coefficients appearing in (4.1) are simply given by whereas if the 

angle dependent redshift is also included then they are given by 
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